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ABSTRACT 

We study explosion characteristics of ultra-stripped supernovae (SNe), which are can¬ 
didates of SNe generating binary neutron stars (NSs). As a first step, we perform stellar 
evolutionary simulations of bare carbon-oxygen cores of mass from 1.45 to 2.0 M 0 un¬ 
til the iron cores become unstable and start collapsing. We then perform axisymmetric 
hydrodynamics simulations with spectral neutrino transport using these stellar evo¬ 
lution outcomes as initial conditions. All models exhibit successful explosions driven 
by neutrino heating. The diagnostic explosion energy, ejecta mass, Ni mass, and NS 
mass are typically ~ 10 50 erg, ~ 0.1 Mq, ~ O.OIMq, and ps 1.3Mq, which are compat¬ 
ible with observations of rapidly-evolving and luminous transient such as SN 2005ek. 
We also find that the ultra-stripped SN is a candidate for producing the secondary 
low-mass NS in the observed compact binary NSs like PSR J0737-3039. 

Key words: binaries: close — stars: evolution — stars: massive — stars: neutron 
supernovae: general — supernovae: individual (SN 2005ek) 


1 INTRODUCTION 

Mergers of binary compact objects, i.e. neutron stars 
(NSs) and black holes (BHs), are promising candidates of 
strong gravitational wave (GW) sources. Event rates of 
these mergers are estimated based on the number of ob¬ 
served binary NSs in our galaxy and population synthe¬ 
sis calculations (e.g., Abadie et al. 2010). These estimates, 
however, have large uncertainty with, roughly speaking, 
two orders of magnitude. Recalling that the compact ob¬ 
jects are formed through gravitational collapse and sub¬ 
sequent supernova (SN) explosions, there should be tran¬ 
sient events generating binary compact objects observ¬ 
able by electromagnetic waves. SN surveys by currently 
working facilities e.g., The Subaru Hyper Suprime-Cam 
(HSC; Miyazaki et al. 2012; Tominaga et al. 2014), Palo- 
mar Transient Factory (PTF; Rau et al. 2009; Law et al. 
2009), Catalina Real-Time Sky Survey (CRTS; Drake et al. 
2009), Panoramic Survey Telescope & Rapid Response Sys¬ 
tem (Pan-STARRSl; Kaiser et al. 2010), and SkyMapper 
(Keller et al. 2007), and also by coming future projects (e.g., 
Large Synoptic Survey Telescope 1 ; LSST) will be able to 

* E-mail: suwa@yukawa.kyoto-u.ac.jp 
1 http://www.lsst.org/lsst/ 


give constraints on the formation rate of transient objects 
including binary compact objects. 

One of the possible candidates for a SN forming a 
close binary system is ultra-stripped SN (Tauris et al. 2015), 
which would launch peculiar type Ib/c SN with a faint and 
fast decaying light curve. Peak luminosity of Type Ib/c 
SNe is determined primarily by the ejected mass of 56 Ni, 
M56 Ni , while the timescale around the peak is determined 
by the diffusion timescale r c oc where M e j is 

the ejecta mass and Ek is the kinetic energy of the ejecta 
(Arnett 1982). Therefore, the low peak luminosity and short 
characteristic time imply the small masses of the ejecta 
and 56 Ni. For instance, type Ic SN 2005ek is one of these 
SNe (Drout et al. 2013; Tauris et al. 2013), whose estimated 
ejecta mass, ~ O(0.1)Mq is notably smaller than typical 
SN Ic, 0(1)M@ (Drout et al. 2011), as well as smaller 56 Ni 
mass, and the explosion energy is also smaller by an order 
of magnitude (O(10 50 ) erg) than typical core-collapse SNe 
(O(10 51 ) erg). To model these rapidly-evolving SNe with 
small ejecta mass, the progenitor stars are thought to be 
stripped much more than canonical stripped-envelope type 
Ib/c SNe, that is, ultra-stripped SNe coined by Tauris et al. 
(2013); Tauris et al. (2015). Besides SN 2005ek and other 
known SNe, ten rapidly-evolving transients were recently 
detected by Pan-STARRSl, which exhibit shorter decay- 
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ing timescale (~ 10 days) than canonical SNe with peak 
barometric luminosities ranging from ~ 10 42 to 10 43 erg 
s -1 (Drout et al. 2014). These ultra-stripped SNe are con¬ 
jectured as products of close binary systems that expe¬ 
rienced strong binary interactions, e.g. common envelope 
phase, which result in close binary compact objects. There¬ 
fore, by assuming that the order of the ultra-stripped SN 
rate is the same as that of the merger rate, we will be able 
to measure the merger rate through SN surveys. The cur¬ 
rent predictions of NS merger rate is between ~ 10 -6 and 
4 x 10 -3 galaxy -1 year -1 (Abadie et al. 2010), which are 
corresponding to ~ 0.01% and 40% of core-collapse SN rate 
(~ 0.01 galaxy -1 year -1 ), while observed ultra-stripped SN 
rate would be ~ 0.1-1% of supernovae (Tauris et al. 2015). 
After the detection of GWs, this conjecture can be tested 
using statistics of GW sources. 

To predict the properties of these transient objects 
for coming LSST era, we investigate the explosion charac¬ 
teristics of the ultra-stripped SNe by means of numerical 
simulations. The current standard model of the explosion 
mechanism for core-collapse SNe is based on a neutrino- 
driven delayed explosion scenario (Colgate & White 1966; 
Bethe & Wilson 1985), in which copious amount of neu¬ 
trinos emitted in the vicinity of a newly-born NS are ab¬ 
sorbed by the surrounding materials and they effectively 
act as heating (see Janka 2012; Kotake et al. 2012; Burrows 
2013; Foglizzo et al. 2015, for latest reviews and references 
therein). In this study, we perform i) stellar evolutionary 
simulations of bare carbon and oxygen (CO) cores with¬ 
out massive hydrogen and helium envelopes, which would 
be possible consequences of close binary interactions such 
as common envelope phase, until iron cores form and ii) ax- 
isymmetric neutrino-radiation hydrodynamics simulations 
to investigate their explosions. The paper is organized as 
follows. Section 2 describes our stellar evolutionary simu¬ 
lations and stellar structure just prior to the collapse. The 
numerical method of following hydrodynamics simulations 
and the results are presented in Section 3. We summarize 
our results and discuss their implications in Section 4. 


2 STELLAR EVOLUTION AND PROGENITOR 
STRUCTURES 

Here, we describe the evolution of CO cores obtained in this 
study. To make initial conditions for hydrodynamics simu¬ 
lations, we first perform stellar evolutionary simulations of 
CO cores with masses of 1.45, 1.5, 1.6, 1.8, and 2.0 Mg sup¬ 
posing that stellar mass loss has already occurred by their 
hypothetical companion NSs. By removing stellar envelope, 
the stellar evolutionary simulations are done with a code 
described in Umeda et al. (2012); Takahashi et al. (2013); 
Yoshida et al. (2014). The nuclear reaction network consists 
of 300 species of nuclei (Takahashi et al. 2013; Yoshida et al. 
2014). Schwarzschild criterion is employed as the convection 
criterion and the convective mixing of the chemical compo¬ 
sition is evolved using diffusion equation. We also take into 
account thermohaline mixing (e.g. Siess 2009) and the diffu¬ 
sion coefficient is adopted from Eq. (2) in Siess (2009). The 
initial chemical compositions of the CO cores are evaluated 
by those after the He burning with the metallicity Z = 0.02. 
Since the C/O ratio of the core depends on the stellar mass, 



Figure 1. Evolutionary path in the plane of the central den¬ 
sity and temperature of C0145 (red line), C016 (green line), 
and CO20 (orange line) models. Labels a and a', b and 6', and 
c,c', and c" correspond to the off-center Ne, O, and Si burnings, 
respectively. Labels a" and b" denote the central Ne and O burn¬ 
ings. 

we assume the mass fractions of 12 C and 16 O, which are 
listed in Table 1. 

It is quite difficult to relate between ZAMS mass and 
CO-core mass in binary systems because mass transfer pro¬ 
ceeds very complicatedly during the evolution. On He star- 
NS binaries, He stars with Mu e < 3.5 Mg and with the ini¬ 
tial orbital period of P or b < 0.5d have an evolutional path 
to ultra-stripped SNe (Tauris et al. 2015). These He stars 
evolved to CO cores with Mco < 1.8 Mg surrounded by very 
thin He envelope. Thus, the mass range of our CO core mod¬ 
els is adequate for the progenitors of ultra-stripped SNe. For 
comparison, we list the stellar mass at the zero-age main- 
sequence (ZAMS), Mzams, which makes the similar CO core 
as shown in Table 1. The mass range of the He cores is 2.6 - 
3.3 Mg. These values are evaluated based on the evolution 
of 9 - 15 Mg stars up to the Ne ignition or the central C ex¬ 
haustion (see also Woosley & Heger 2015). We assume the 
convective overshooting as a diffusion process until the core 
He burning in these calculations. The diffusion coefficient is 
adopted from Eq. (2) in Herwig (2000) and the parameter 
on the scale height / is set to be 0.01. 

We continue evolutionary simulations until the iron 
cores become unstable and start collapsing. Here we briefly 
explain the evolution of our models. More details will be 
shown in Yoshida et al. (2015, in preparation). Figure 1 
shows the evolution of the central temperature ( Tc ) and 
density (pc) of C0145, C016, and CO20 models. Here, the 
units of Tc and pc are Kelvin and g cm -3 , respectively. 
Figure 2 is the Kippenhahn diagram of C0145 model, which 
displays the time evolution of convective regions and the en¬ 
ergy generation inside the star. Although all models finally 
form Fe cores, the evolution paths depend on the CO-core 
mass. In C0145 model, off-center Ne burning is ignited (Ne- 
flash) after several C-shell burnings (label a in Figures 1 and 
2). The following O-burning is also ignited at an off-center 
region (label b ). Then, the burning front gradually moves 
into the center and the temperature at the front increases. 
When the central density reaches log pc ~ 8.2, Si-burning 
is ignited at the burning front and the convective layer ex¬ 
tends to M r ~ 0.9Mq (label c). The burning front continues 
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Time to Collapse (year) 

Figure 2. Kippenhahn diagram of 00145 model. Green-hatched 
areas indicate convective regions. Labels a, b, and c correspond 
to the off-center Ne, O, and Si-burnings at 17.1 years, 10.7 years, 
and 38.8 days before core-collapse. 


to move inwards and the O/Ne-enriched material in the cen¬ 
tral region is burned to Fe-peak elements. When the central 
density reaches log pc ~ 8.5, the burning front reaches the 
center and the central temperature raises steeply. The Fe 
core grows up to ~ 0.5Mg after the burning front reaches 
the center. Then, shell Si-burning occurs and the Fe core 
grows up further, and finally, the core collapses. 

CO 16 model evolves on a different evolution track as 
shown in Figure 1. After the Ne flash (label a') and off- 
center O burning (label b’), the burning front reaches the 
center. As a result, the Si core with 0.6 Mq is formed and it 
grows up by O-shell burnings. Then, the Si flash occurs at 
M r = 0.03 Mq (label c 1 ). The Si flash forms large convective 
region up to M r ~ 1.1 Mq above the burning front. After 
that, Si-shell burnings also follow. The burning front again 
gradually moves inwards and finally reaches the center after 
the second Si-shell burning. Finally, the Fe core grows to 1.34 
Mq and collapses. The evolution of C015 model is similar 
to that of C016. The Ne-flash and Si-flash occur at higher 
central density. 

The evolution of CO20 model is similar to those of mas¬ 
sive stars that collapse to normal SNe until O-shcll burning. 
Ne and O burnings are ignited at the center and the con¬ 
vective core forms (labels a" and b"). The Ne and O-shell 
burnings extend the central Si core after O-core burning. 
However, Si burning starts at an off-center region when the 
central density reaches log pc ~ 8.6 (label c"). The off-center 
Si burning forms a large convective region to 1.2 Mq. The 
burning front moves inwards and reaches the center during 
Si-shell burning. Finally, an 1.37 Mq Fe core is formed, and 
then, it collapses. The evolution of C018 model is similar 
to CO20 model, i.e., Ne and O are ignited at the center and 
Si ignites at an off-center region. The central density at the 
ignition of off-center Si-burning is higher than that of CO20. 

We note that recently the evolution of 9 - 11 Mq stars 
has been investigated in Woosley & Heger (2015). Since they 
considered artificial energy deposition to the cooler under¬ 
lying zone of the off-center burning front, the burning front 
rapidly moved inward. We do not consider this effect, but we 
instead set larger coefficient value of thermohaline mixing. 
In CO 145 model, the central contraction makes the temper¬ 


ature at the burning front high enough to ignite Si before the 
front reaches the center. Even so, CO 145 model forms an Fe 
core. Although the inward motion of the burning front de¬ 
pends on the instability at the base of the burning front and 
has uncertainty, we expect that the stars igniting Ne-flash 
will form an Fe core and collapse. 

The mass fraction distributions of C0145 and C016 
models at the final step of the stellar evolution simulations 
are shown in Figure 3. These models consist of a 1.33 and 
1.34 Mq Fe cores surrounded by thin Si, O/Si, O/Ne, and 
O/C layers. Composition-inverted layers have been seen dur¬ 
ing the Ne-flash and Si flash. The composition inversion has 
been removed during the inward motion of the burning front. 
We do not see large differences in Fe core masses in our CO 
core models (see Table 1). 

We assume the mass fractions of C and O in the CO 
cores based on the results of single star models. In the case 
of binary evolution, ultra-stripped CO cores have lost the H- 
rich envelope before or during the early phase of the He-core 
burning. Stripping the H-rich envelope prevents the growth 
of He core and shortens the time scale of the He-core burn¬ 
ing. Thus, the C/O ratio of an ultra-stripped CO core is 
larger than the CO core in single stars (Brown et al. 2001). 
Wellstein & Langer (1999) investigated the evolution of var¬ 
ious binary systems of massive stars. They obtained that a 
binary system consisting of 13 Mq and 12 Mq stars with 
an initial orbital period of 3.1 days experiences Case B+BB 
mass transfer and the primary star becomes a 1.42 Mq He 
star with CO core mass being 1.31 Mq and the central car¬ 
bon mass fraction being 0.40. 

We calculate the evolution of 1.45 and 1.6 Mq CO core 
models with a large C/O ratio to investigate the C/O-ratio 
dependence of ultra-stripped SN progenitors. We set the 
mass fractions of C and O in these models to be 0.471 and 
0.500, respectively. Properties of the models are shown as 
CO 145c and CO 16c in Table 1. We do not see large system¬ 
atic differences between normal and large C/O-ratio models. 
CO 145c model has an Fe core slightly smaller than CO 145, 
while the Fe core of CO 16c is slightly larger than that of 
CO 16 for the criterion that the Fe core is determined by the 
electron fraction smaller than 0.495. The difference of the 
compactness parameter between the normal model and large 
C/O-ratio model is within 10 %. We consider that properties 
of ultra-stripped SN progenitor do not strongly depend on 
the C/O ratio. It should be noted that the evolution of the 
central region of the CO cores depends on the C/O ratio. 
The convective regions of the central or shell C-burnings be¬ 
come large in the models of the large C/O ratio. In CO 145c 
model, Ne is ignited at 0.35 Mq in the mass coordinates. The 
Ne/O-burning front reaches the center before the Si ignition. 
Then, Si is ignited at 0.04 Mq. In C016c model, Si-burning 
occurs at the center as a flash and, then, it turns to steady 
burning. Nevertheless, we do not see large difference by the 
C/O ratio in the region where the structures and compo¬ 
sition are mainly determined by the Si-shell burning. Once 
the Fe core is formed by the central or off-center Si-burning, 
the core grows through the following Si-shell burning. Thus, 
we consider that properties of SN progenitors such as Fe 
core mass and compactness parameter are insensitive to the 
C/O ratio. 

These CO cores have a thin CO envelope where shell 
C burnings scarcely affect the composition. We listed the 
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Figure 3. Mass fraction distributions of C0145 (left) and 0016 (right) models. Red, black, sky-blue, blue, pink, yellow-green, green, 
and orange lines correspond to 1 H, 4 He, 12 C, ls O, 20 Ne, 24 Mg, “Si”, and “Fe”. “Si” and “Fe” indicate the element groups in Z = 14 — 21 
(Si-Sc) and Z > 22 (Ti-), respectively. The scale of the horizontal axis changes at 1.3 Mg. 



Figure 4. Density structure of CO cores at precollapse phase as a 
function of mass coordinate. Red, blue, green, purple, and orange 
lines represent models with CO cores of 1.45, 1.5, 1.6, 1.8, and 
2.0M 0 , respectively. 

masses of the envelope, M env , which the convective region 
did not reach during C-shell burnings, and their binding en¬ 
ergy, -Ebind.env, in Table 1. The envelope mass is less than 
0.01 M 0 for all the models. The binding energy of the enve¬ 
lope is (0.4 — 2) x 10 48 erg, which is smaller than one per cent 
of the binding energy of the whole star (f?bind, see Table 1). 

Although we do not take into account He-rich envelope, 
the progenitors of ultra-stripped SNe in NS binary systems 
could have the envelope (Tauris et al. 2013; Tauris et al. 
2015). Tauris et al. (2015) showed the dependence of the 
He mass and the binding energy of the He-envelope on the 
final CO core mass. Here, we estimate the binding energy of 
the He-envelope from Figure 15 in Tauris et al. (2015) and 
from their discussion that the He mass is smaller than 0.2 
M 0 (see §4.1.1 in Tauris et al. 2015). The estimated range 
of the envelope is listed in Table 1. The binding energy of 
the He envelope is at most several times 10 49 erg. Thus, we 
expect that the envelope has an small contribution in the 
binding energy for ultra-stripped SNe. 

Figure 4 presents density structures at the end of evo¬ 
lutionary simulations. More detailed properties are summa¬ 
rized in Table 1. Following O’Connor & Ott (2011), we de¬ 


fine the compactness parameter, £m, by 

M/Mg 

? _R(M)/1000 km’ 1 j 

where R(M ) is the radius with M being mass enclosed in¬ 
side R. We employ £ 1 . 4 , which is measured by M r = 1.4Mg. 
Note that other groups have employed different values such 
as £ 2 . 5 , £ 1 . 75 , £ 1.5 (O’Connor & Ott 2011; Ugliano et al. 
2012; O’Connor & Ott 2013; Sukhbold & Woosley 2014; 
Nakamura et al. 2014; Pejcha & Thompson 2015). Since our 
models have small CO core mass, we rely on a differ¬ 
ent parameter. Note also that in the original definition 
by O’Connor & Ott (2011), £ is measured at the time of 
bounce. Later, Sukhbold & Woosley (2014) suggested that 
£ measured just prior to collapse correlates with £ at bounce. 
Hence, in this paper we evaluate £ with the precollapse den¬ 
sity structures. We show two different definitions of the iron 
core mass, Mf 6 : one is based on the electron fraction and the 
other on mass fraction of heavy elements. The iron core mass 
does not depend strongly on CO core mass since, broadly 
speaking, the core collapse sets in when the iron core mass 
exceeds Chandrasekhar mass (Baron & Cooperstein 1990). 

3 HYDRODYNAMICS SIMULATIONS 
3.1 Numerical methods 

For the hydrodynamics simulations, we employ a two- 
dimensional (2D) neutrino-radiation hydrodynamics code, 
which is developed and used for investigating SN explo¬ 
sion mechanism (Suwa et al. 2010, 2011, 2013; Suwa et al. 
2014; Suwa 2014). With ZEUS-2D code (Stone & Norman 
1992) as a base for the solver of hydrodynamics, we em¬ 
ploy an equation of state of Lattimer & Swesty (1991) with 
an incompressibility K = 220 MeV and solve the spectral 
transfer of neutrinos by the isotropic diffusion source ap¬ 
proximation (IDSA) scheme (Liebendorfer et al. 2009) that 
splits the neutrino distribution function into two compo¬ 
nents, both of which are solved by using separate numerical 
techniques. Weak interaction rates are implemented follow¬ 
ing Bruenn (1985). We solve transfer of electron-type neu¬ 
trinos (u e ) and antineutrinos (p e ), but heavier leptonic neu- 
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Table 1. Properties of progenitor models 


Model 

Mqo 

[Mq\ 

AfzAMS ° 
[M®\ 

N C 

No 

Radius 

[10 4 km] 

$1.4 

Mf, 6 
[Mq\ 

Mf« c 
[M©] 

Me nv 

[io- 3 m 0 ] 

^bind,env 

[10 48 erg] 

jpd 

^Dind,He 
10 49 erg] 

-^bind 
[10 50 erg] 

C0145 

1.45 

9.75 

0.360 

0.611 

1.28 

0.468 

1.33 

1.32 

4.80 

2.28 

0.1-1.3 

3.97 

C0145c 

1.45 

9.75 

0.471 

0.500 

1.78 

0.502 

1.31 

1.28 

1.45 

0.46 

0.1 1.3 

3.99 

0015 

1.5 

10.0 

0.360 

0.611 

2.20 

0.600 

1.31 

1.29 

1.60 

0.38 

0.5-1.7 

4.17 

C016 

1.6 

10.5 

0.350 

0.621 

1.93 

0.912 

1.34 

1.34 

2.63 

0.80 

0.7-2.8 

5.39 

C016c 

1.6 

10.5 

0.471 

0.500 

1.73 

0.850 

1.36 

1.34 

2.42 

0.95 

0.7-2.8 

4.69 

0018 

1.8 

11.5 

0.350 

0.621 

2.64 

0.851 

1.35 

1.37 

8.09 

2.20 

1.2-3.8 

4.92 

CO20 

2.0 

12.8 

0.330 

0.641 

3.66 

0.968 

1.37 

1.37 

8.48 

1.66 

2.2-4.7 

5.72 


“ This is corresponding mass of ZAMS stars, which make the same mass of CO core. 
b The iron core mass is determined by the mass with the electron fraction of Y e < 0.495. 

c The iron core mass is determined by the mass with the mass fraction of the element groups in Z ^ 22 (Ti-) larger than 0.5. 
d These values are taken from Figure 15 in Tauris et al. (2015). 


trinos are not taken into account. In our 2D simulations, 
axial symmetry is assumed and “ray-by-ray-plus” approach 
is implemented for multi-dimensional treatment of neutrino 
transfer (Buras et al. 2006). Spherical coordinates (r, 9) with 
logarithmic zoning in the radial direction and constant zon¬ 
ing in 9 are used. The simulations are performed on a grid 
of 300 radial zones extending up to 5000 km with the small¬ 
est grid width being 1 km at the center and 128 equidistant 
angular zones covering 0 < 9 < n for 2D simulations. For 
neutrino transport, we use 20 logarithmically spaced energy 
bins ranging from 3 to 300 MeV. 

Note that at the current moment there is no complete 
model for the explosion mechanism of core-collapse SNe de¬ 
spite the long-lasting efforts. Recent multi-dimensional sim¬ 
ulations of neutrino radiation hydrodynamics have shown 
diversity of numerical results (e.g. Bruenn et al. 2013; 
Suwa et al. 2014; Dolence et al. 2015; Melson et al. 2015b, 
used the same progenitor models and obtained different re¬ 
sults). In addition, most exploding simulations exhibited or¬ 
der of magnitude smaller explosion energy than observation 
and remnant compact objects often gained mass above the 
maximum mass of a NS, which would collapse to BHs in¬ 
stead of NSs. 

We speculate, however, that all these implementations 
could universally reproduce ultra-stripped SNe. The reason 
for this speculation is that typical explosion energy of the 
ultra-stripped SNe is smaller (O(10 50 ) erg) than canonical 
SNe (O(10 51 ) erg), that is, stellar structures of progenitors 
for these different kinds of SNe are different. The explosion 
energy is determined naively by the binding energy of pro¬ 
genitor layers in the vicinity of mass cut (remnant compact 
object) so that stars with small binding energy are possi¬ 
ble candidates of the weak explosions. This hypothesis was 
already applied for explanation of type Iln-P SNe like SN 
2009kn, which exhibit narrow emission lines, a short plateau 
phase in light curve, and small amount of 56 Ni. These obser¬ 
vational features were reproduced by an electron-capture SN 
model, in which the progenitor star consists of an O-Ne-Mg 
degenerated core and a very thin envelope, giving the small 
explosion energy of ~ 10 50 erg (Moriya et al. 2014 based on 
the explosion simulation by Kitaura et al. 2006). As we will 
show in the following subsection, we obtain explosion en¬ 
ergy of O(10 50 ) erg, which is the same order of magnitude 
as the binding energy exterior to a remnant compact object 



Mass coordinate [M . ] 


Figure 5. Density profiles at 5 ms after bounce. Line colors show 
the same models as Figure 4. 


at precollapse phase. We infer that even if we could per¬ 
form a realistic simulation that would reproduce canonical 
explosion energy for canonical SN progenitors, the results 
obtained with weakly bound progenitors will be a weak ex¬ 
plosion with its energy ~ 10 50 erg. Since the explosion en¬ 
ergy of such simulations, however, could be larger than that 
of ours, it is safe to consider that our results give at least a 
lower limit of the explosion energy. 

At first, we perform spherically symmetric (ID) simu¬ 
lations up to 5 ms after the bounce, which is determined by 
the largest central density, and 2D simulations follow them. 
Note that all ID simulations fail to explode even when we 
continue simulations until several hundred milliseconds after 
the bounce. 2 
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3.2 Results 

We show the density profiles at 5 ms after the bounce in Fig¬ 
ure 5. One can find that the density profiles of M r < I.IMq 
depend only weakly on the CO core mass, but above this 
mass coordinate they differ from each other due to the dif¬ 
ferent progenitor structure. This difference leads to different 
evolution of shocks, which is described in the following. Note 
that the density at a certain mass coordinate above ~ I.IMq 
monotonically increases with the CO core mass. 

Figure 6 presents evolution of entropy (left halves) and 
radial velocity (right halves) distributions at 100 (left top 
panel), 200 (right top), 250 (left bottom), and 350 (right 
bottom) ms after the bounce for model C015. After core 
bounce, convective motion sets in initially due to an unstable 
entropy configuration (prompt convection) and subsequently 
due to neutrino heating in the gain region (neutrino-driven 
convection), but shock morphology is maintained to be ap¬ 
proximately spherically symmetric (see left top panel). As 
the turbulent motion is further enhanced by the neutrino 
heating, the shock is deformed from spherical symmetry (see 
right top). Once the postshock pressure (thermal pressure 
and turbulent pressure) overwhelms the ram pressure above 
the shock, the shock starts expanding and an explosion is 
launched. In this model, the expanding shock is rather spher¬ 
ical since the small envelope mass results in a small mass 
accretion rate onto the shock and small-scale convection 
dominates over large-scale motion driven by the standing 
accretion shock instability (SASI; Blondin et al. 2003). This 
leads to a small kick velocity of protoneutron stars (PNS). 
One can observe a cold downflow onto the PNS even after 
the onset of the explosion (see an isolated blue region in the 
bottom right panel of Figure 6), which could increase the 
PNS mass. This downflow, however, has a small solid angle 
and the mass accretion rate is considerably small. Thus, the 
mass accreting onto the PNS at this moment is negligible. 
The PNS mass evolution will be discussed later. 

Shock evolution processes in all 2D simulations are pre¬ 
sented in Figure 7. One can see that the three small mass 
models (C0145, C015, and 0016) explode approximately 
at identical time, while more massive cores show later explo¬ 
sion. This is a consequence of the different envelope mass on 
the iron core. The early exploding models have clear phase 
transition from a slowly expanding phase to a rapidly ex¬ 
ploding phase at ~ 200 ms after the bounce. On the other 
hand, late exploding models show oscillations of the shock 
radius before the explosion. This is due to convection and 
SASI. The difference in the onset time of the explosion is a 
result of different mass accretion rate evolution due to the 
different envelope structure (see Figure 5). The later onset 
of the explosion leads to larger PNS mass as shown in Figure 
8. Here we define PNS as the region with density above 10 11 
g cnF 3 . C0145 and C015 models form a PNS of baryonic 
mass « 1.35Mq, while other models give larger PNS mass. 

In Table 2, we summarize results of our hydrodynamics 

2 Note that recent studies of full three-dimensional (3D) hydro¬ 
dynamics simulations showed that 2D simulations are more favor¬ 
able for explosion than those of 3D (see e.g., Hanke et al. 2013; 
Takiwaki et al. 2014) so that we should note that the 2D simula¬ 
tions would optimize the conditions for successful explosions (but 
see also Melson et al. 2015a). 


simulations. tg na i denotes the final postbounce time of each 
simulation. The quantities listed in other columns are all 
measured at ffi na i. -Rsh is the angle-averaged shock radius, 
Eexp is diagnostic explosion energy, which is defined as the 
integral of the sum of specific internal, kinetic and gravita¬ 
tional energies over all zones with positive value, MNS.baryon 
is baryonic mass of the remnant NS calculated by integration 
over grid of p > 10 11 g cm -3 , MNS.grav is the corresponding 
gravitational mass, M e j = Mco — MNS.baryon is the ejecta 
mass, M56 Ni is mass of 56 Ni, and fkick is the estimated kick 
velocity of NSs. Note that these quantities are not the final 
outcome of the simulations, since all the simulations were 
terminated before the system relaxes to a stationary state 
to save the computational time. The gravitational mass is 
calculated by the baryonic mass using the following relation 
(Timmes et al. 1996) 

Mbaryon _ Mg ra v _ (Mgrav \ 

Mq Mq \ Mq ) 

56 Ni mass is calculated using tracer particle method (e.g., 
Nagataki et al. 1997). We assume that the mass elements 
with the maximum temperature being over 5 x 10 9 K achieve 
nuclear statistical equilibrium and synthesize 56 Ni com¬ 
pletely. This gives just an approximate estimate. For more 
realistic calculation we need to perform detailed nucleosyn¬ 
thesis calculation, which is beyond the scope of this paper. 
The NS kick velocity is estimated by assuming the linear 
momentum conservation of the whole progenitor star, i.e., 
assuming that anisotropic mass ejection leads to NS kick 
(e.g., Wongwathanarat et al. 2013). The linear moment of 
ejecta is calculated by 

Pej = [ P V dV, (3) 

J p< 10 11 g cm -3 ,i) r >0 

where v is the velocity vector and v r is its radial com¬ 
ponent. The NS kick velocity is then given by Vkick = 
—P e j/MNS.baryon- Since the axial symmetry is assumed in 
our simulations, the kick velocity may be overestimated due 
to the existence of preferable direction of NS kick, i.e. sym¬ 
metry axis. Additionally, the stochastic nature of postshock 
turbulent flow would also change the degree of asymmetry 
of ejecta so that the the initial small perturbation could 
change the kick velocity significantly (Scheck et al. 2006). 
More statistical study is needed to pin down this issue. It 
can be argued that small envelope, not small iron core itself, 
which can rapidly accelerate shock, would generally lead to 
small kick velocity due to too short time for SASI to build up 
(see also, e.g., Podsiadlowski et al. 2004; Bogomazov et al. 
2007). 


4 SUMMARY AND DISCUSSION 

We have performed both stellar evolution simulations of bare 
CO cores and explosion simulations for the end product of 
the CO cores for modeling ultra-stripped type Ic SNe. We 
have found that all CO cores with mass from 1.45 to 2.0 
Mq resulted in explosion with energy of O(10 5 °) erg, which 
left NSs with gravitational mass from ~ 1.24 to 1.44 Mq 
and ejecta from ~ 0.1 to 0.4 Mq with synthesized 56 Ni 
of O(10 _2 )Mq. These values are compatible with obser¬ 
vations of ultra-stripped SN candidates (Drout et al. 2013; 
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Figure 6. Specific entropy (in unit of fcg baryon -1 ; left halves of the individual panels) and radial velocity (in unit of 10 4 km s 1 ; right 
halves) profiles at 100 (left top panel), 200 (right top), 250 (left bottom), and 350 (right bottom) ms after the bounce for model C015. 
In the entropy plots, blueish (reddish) colors represent small (large) entropy. In velocity plots, red region is expanding (positive radial 
velocity) and blue region is accreting (negative radial velocity). 


Tauris et al. 2013; Tauris et al. 2015). For SN 2005ek, M e j « 
0.2-0.7 Mg and Mni « 0.02-0.05 M 0 are appropriate to fit 
its light curve. The event rate of these SNe is estimated 
as ~ 1% of core-collapse SN rate (Drout et al. 2013, 2014), 
which is also compatible with an NS merger rate estimation 
(Abadie et al. 2010). 

We took a different approach from previous studies on 
ultra-stripped SNe (Tauris et al. 2013; Tauris et al. 2015). 
In previous works, they self-consistently performed stellar 
evolutionary simulations until oxygen burning phase with 
self-consistent mass loss driven by wind but explosion calcu¬ 
lations were based on phenomenological modeling with three 
free parameters; kinetic energy of SN, Ni mass, and mass cut 
(i.e. NS mass). Based on this model, they found that ultra- 
stripped SN model could account for the light curve of SN 


2005ek quite well. In our work, on the other hand, we per¬ 
formed stellar evolutionary simulations until the last phase 
of evolution, i.e., iron core collapse, but for initially bare 
CO cores without mass loss. For the explosion phase, we 
performed neutrino radiation hydrodynamics simulations to 
calculate explosion energy, Ni mass, and NS baryon mass in 
self-consistent manner. In this sense, this work is comple¬ 
mentary to previous works. We found that outcomes of our 
hydrodynamics simulations for C0145, C015, and C016 are 
close to model parameters of Tauris et al. (2013) that well 
fit light curve of SN 2005ek (Mns = 1.3 Mq, M e j = 0.2M 0 , 
E exp = 5 x 10 50 erg, and Mm = O.O5M 0 ). The light curve 
calculated with our hydrodynamical models will be pre¬ 
sented in the forthcoming paper. 

For C0145, CO!5, CO!6, and CO!8 models, the baryon 
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Table 2. Summary of simulation results 


Model 

^final 

[ms] 

R s h b 

[km] 

Eexp C 

[B] 

-^NS, baryon 

[Mq] 

-^NSjgrav 6 

[Mq] 

M ej 7 
[10 -1 M e ] 

M m 9 
[10 -2 M q ] 

^kick 
[km s -1 ] 

00145 

491 

4220 

0.177 

1.35 

1.24 

0.973 

3.54 

3.20 

0015 

584 

4640 

0.153 

1.36 

1.24 

1.36 

3.39 

75.1 

C016 

578 

3430 

0.124 

1.42 

1.29 

1.76 

2.90 

47.6 

0018 

784 

2230 

0.120 

1.49 

1.35 

3.07 

2.56 

36.7 

0020* 

959 

1050 

0.0524 

1.60 

1.44 

3.95 

0.782 

10.5 


a The final time of simulations measured by postbounce time. 
b The angle-averaged shock radius at ffl na i. 

c The explosion energy in unit of B (=10 51 erg) at £fi na i, which is still increasing. 

The baryonic mass of NS at £fi na i. 
e The gravitational mass of NS computed by Eq. (2) at ifl na i. 
f The ejecta mass at tfi na i. 

9 The Ni mass at t(\ rl a i. 
h The kick velocity at ifj na j. 

1 Note that this model is marginally exploding. 



Time after bounce [ms] 

Figure 7. Evolutions of the radius of shocks. 



Time after bounce [ms] 

Figure 8. Time evolutions of PNS mass (define by p > 10 11 g 
cm -3 ). 

mass of the remnant NSs is in the range between 1.35Mq 
and 1.49Mq. The corresponding gravitational mass are 1.24 
- 1.35 Mq for this baryon mass range depending only weakly 
on equations of state for neutron stars hypothetically em¬ 
ployed. These values for the gravitational mass agree well 


with those for secondary NSs in observed compact binary 
systems (Lattimer 2012). This suggests that such secondary 
NSs may be formed from ultra-stripped SNe. 

We showed that 1.45 — 2.0 Mq CO star models ignite 
Ne at the center or an off-center region. They form an 
Fe core and none of them evolve to electron-capture SN. 
Thus, the upper-limit of the CO core mass for electron- 
capture SN should be less than 1.45 Mq. Tauris et al. (2015) 
adopted MoNe.f = 1.43 Mq as an approximate upper limit 
for electron-capture SN. In the case of single star evolution 
in Takahashi et al. (2013), a 10.8 Mq model ends its evolu¬ 
tion as an electron-capture SN and an 11.0 Mq model ignites 
Ne at an off-center region. Off-center Ne ignition and grad¬ 
ual increase in the central temperature around the central 
density pc ~ 10 9 g cm -3 could be predictions of the Fe core 
formation and core-collapse SN. The evolution of CO cores 
less massive than 1.45 Mq will be shown in the forthcoming 
paper (Yoshida et al. 2015, in preparation). 

Finally, we discuss the eccentricity of the binary system 
formed after the SN that leaves a secondary NS. Due to the 
mass ejection from the system, the binary system obtains 
the eccentricity, e, after the SN explosion. The eccentricity 
can be evaluated by (Bhattacharya & van den Heuvel 1991) 


Mi - M{ 

6 ~ M{ + Mns ’ 


where Mi is the mass of exploding star before (after) the ex¬ 
plosion indicated as i (/) and Mns is the mass of the primary 
NS. Here, we assumed that the mass ejection occurs quickly 
and during the explosion the positions of these stars do not 
change. By giving the ejected mass, M e j = Mj — M{ = 
0.3Mq and M( = Mns = 1.3 Mq, we get e ~ 0.12, which 
is compatible with one of observed binary NSs, J0737-3039, 
whose current eccentricity is 0.088 and estimated eccentric¬ 
ity at birth of second pulsar is 0.11 (Piran & Shaviv 2005). 
The small center of mass velocity of this system also implies 
a small ejecta mass and slow pulsar kick (DalPOsso et al. 
2014). 
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